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Summary. — The relative-space-time-transformation (RSTT) paradigm and the 
interpretation of the burst-structure (IBS) paradigm are applied to probe the origin 
of the time variability of GRBs. Again GRB 991216 is used as a prototypical case, 
thanks to the precise data from the CGRO, RXTE and Chandra satellites. It is 
found that with the exception of the relatively inconspicuous but scientifically very 
important signal originating from the initial "proper gamma ray burst" (P-GRB), 
all the other spikes and time variabilities can be explained by the interaction of the 
accelerated-baryonic-matter pulse with inhomogeneities in the interstellar matter. 
This can be demonstrated by using the RSTT paradigm as well as the IBS paradigm, 
to trace a typical spike observed in arrival time back to the corresponding one in the 
laboratory time. Using these paradigms, the identification of the physical nature of 
the time variablity of the GRBs can be made most convincingly. It is made explicit 
the dependence of a) the intensities of the afterglow, b) the spikes amplitude and c) 
the actual time structure on the Lorentz gamma factor of the accelerated-baryonic- 
matter pulse. In principle it is possible to read ofi' from the spike structure the 
detailed density contrast of the interstellar medium in the host galaxy, even at very 
high redshift. 



It is well known that one of the most successful cognitive tools in relativistic astro- 
physics has been the analysis of the time structure of signals received at a variety of 
wavelengths. Time variabilities, however, have not always been of significance in rela- 
tivistic astrophysics. In the case of pulsars, for example, only the period of the average 
pulsar signal and its monotonic lengthening with time have been essential in identify- 
ing pulsars as rotating neutron stars (Hewish, Bell, et al. (1977) [9]). Furthermore, the 
modulations of the pulsar signal, periodic in time, have been essential for the identifi- 
cation of binary pulsars and to give the first evidence for gravitational waves (Hulse & 
Taylor (1975) [11]). The secular variation of the shape of the pulse yields information 
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relating the role of the spin of the neutron star and its angular momentum to additional 
general relativistic effects (Damour & Ruffini (1974) [7], Kramer (2001) [12]). Even in 
this very successful example, there is a broad range of effects connected to pulsars whose 
role in relativistic astrophysics and fundamental physics is null. We quote, as an exam- 
ple, the well known time delay in the arrival time of pulsar signals, inversely proportional 
to the square of the radiation frequency of observation, sec e.g. Rccs. ct al. (1975) [17] 
page 26. Such an effect is not due to the mass of the photon, as one might have hoped 
for fundamental physics reasons (see e.g. Ohanian & Ruffini (1994) [10], page 117), but 
simply due to the dispersion by the electrons in interstellar plasma. 

In approaching the analysis of GRB signals it is similarly essential to untangle infor- 
mation about the astrophysical system producing the GRBs, which is certainly in the 
realm of relativistic astrophysics, from other parts of the signal, also of similar magni- 
tude and structure, which can instead be traced back to the environment in which the 
astrophysical process occurs and in that sense may very well belong to the domain of 
classical astronomy. 

In the last decade an enormous number of papers have been written trying to link 
all the structures observed by the BATSE experiment on the CGRO satellite to the 
intrinsic properties of an unknown GRB source, whose properties should be determined, 
hopefully, by those observations. In the "internal shock models" of GRBs, which are 
currently very popular, it is assumed that every spike in the burst in the range At ~ 1 
s to ~ 50 s is directly related to the physical properties of the "inner engine" (sec; e.g. 
Piran (2001) [13] and references therein). The fact that it is difficult to explain the long 
bursts has led the theorists working on the "internal shock model" to introduce a new 
family of models, in which the source of GRBs has a prolongued action in time. We shall 
see below that a simpler and very different explanation can be found. 

The electromagnetic black hole (EMBH) model (Ruffini (1998) [18], Preparata, et 
al. (1998a) [15], Preparata, et al. (1998b) [16]) relates the origin of the energy of GRBs 
to the extractable electromagnetic energy of an EMBH [5] via the vaciiTim polarization 
process occurring during the gravitational collapse leading to the formation of an EMBH 
[8]. The first step in this theory is the definition of the dyadosphere [18, 15], an extended 
region outside the EMBH horizon formed of an optically thick plasma of electron-positron 
pairs and radiation whose energy E^ya is related to the mass fi = M/Mq and electro- 
magnetic parameter ^ = Q/ (mVg) of the EMBH by the relation: 



(1) 



Ed: 



ya 



2r+ 



1 - ^ 



1 - ^ 



where r+ = 1.47 x 10^/Lt(l -|- yT— ^) is the horizon radius and rds = 1-12 x lO^vT*? is 
the dyadosphere radius and, as usual, M and Q are the mass-energy and charge of the 
EMBH and G is the Newton constant of gravity. 

The evolution of this pair-electromagnetic plasma leads to the formation of a sharp 
pulse (the PEM pulse) that very rapidly reaches a Lorentz gamma factor of 10^ and 
higher. The subsequent interaction of this pulse with the baryonic matter of the remnant, 
left over from the gravitational collapse of the protostar, and with the interstellar medium 
(ISM) leads to the different eras of the GRBs. It is useful to parametrize the baryonic 
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Fig. 1. - The Lorentz gamma factor corresponding to the different eras of GRB991216 is given 
as a function of the radial distance from the EMBH. Details in Letter 1. 



mass Mb of the remnant by introducing the dimensionless parameter B: 
(2) Mbc^ = BEdya. 

The confrontation of the theoretical model with the observational data allows an 
estimate for the values of the EMBH parameters. It also allows us to probe the density 
of the baryonic material in the remnant, in the ISM as well as in the stellar distribution 
within a few parsecs of the EMBH (see [22, 23]). 

In RufBni, ct al. (2001a) [21] we presented the rolative-space-time-transformation 
(RSTT) paradigm, leading to the diagram relating the Lorentz gamma factor of the 
pulse to the space and time parametrization both in the comoving and in the labora- 
tory frame for the case of GRB 991216 (sec Fig. 1). In Ruffini, ct al. (2001b) [22] wc 
introduced the interpretation of the burst-structure (IBS) paradigm, presenting a dras- 
tic separation between the proper-gamma-ray burst (P-GRB) and the E-APE, the "not 
burst component of the GRB", sec Fig. 2. 

It is important to stress that the results obtained in the IBS paradigm are of general 
validity for a variety of GRB sources based on a single gravitational collapse event. What 
makes the EMBH model uniqueness testable are: 

• the energetics [18], 

• the time structure of the P-GRB [25], 
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Fig. 2. - The flux of the afterglow of GRB991216, as computed using the best fit of the data 
obtained in Letter 2, is given as a function of the arrival time. The dashed (dotted) hne 
corresponds to Eq.(6) (Eq.(9)) in the text. 



• the spectral information of the P-GRB [3] . 

Once again we use GRB 991216 as a prototypical case due to the excellent data from 
the BATSE (BATSE Rapid Burst Response (1999) [1]), RXTE (Corbet & Smith (2000) 
[6]) and CHANDRA (Piro, et al. (2000) [14]) satellites, although the conclusions will be 
applicable to all GRBs. In a certain sense, in this paradigm all the features of a GRB 
are divided in two very distinct phases: 

• the first, prior to decoupling and ending with emission of the P-GRB, which we 
shall call the "injector phase" ; 

• the second, the "beam-target" phase, in which the accelerated-baryonic-matter 
(ABM) pulse (the beam) interacts with the interstellar medium (the target). 

The P-GRB, which is clearly identifiable in the enlargement in Fig. 3a, is emitted when 
the condition of transparency is reached by the pair-electromagnctic-baryonic matter- 
pulse, the PEMB pulse. As already stressed RufRni, et al. (2001b) [22], the properties 
of the P-GRB are directly linked (Ruffini, et al. (1999) [19], Rufiini, et al. (2000) [20], 
Bianco, et al. (2001) [2], Ruffini, et al. (2001g) [25]) to the internal properties of the GRB 
source and to the detailed structure and energy distribution in the dyadosphere of the 
EMBH (Ruffini (1998) [18], Preparata, et al. (1998a) [15], Preparata, et al. (1998b) [16]). 
These results are essential in identifying the extractable energy of the EMBH, introduced 
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Fig. 3. - a) Flux of GRB 991216 observed by BATSE. The enlargement clearly shows the P-GRB 
(see Letter 2). b) Flux computed in the collision of the ABM pulse with an ISM cloud with 
the density profile given in Fig. 4. The dashed line indicates the emission from an uniform ISM 
with n = lcm~^. The dotted line indicates the BATSE noise level. 
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in Christodoulou & Ruffini (1971) [5], as the GRB energy source. The operational tool of 
the energy extraction process is the vacuum polarization process introduced in Damour 
& RufRni (1975) [8]. Similarly the intensity ratio of the P-GRB to the afterglow gives a 
precise measurement of the matter left in the remnant during the process of gravitational 
collapse of the progenitor star to the EMBH, see Ruffini, et al. (2000) [20] and Ruffini, 
et al. (2001b) [22] and references therein. 

All the above results clearly involve relativistic astrophysics. Let us now turn to the 
afterglow and apply the RSTT and IBS paradigms in order to understand its detailed 
time structure. 

The afterglow is emitted as the ABM pulse plows through the interstellar matter 
engulfing new baryonic material (RufRni, et al. (20011) [24]). In our previous works we 

were interested in explaining the overall energetics of the GRB phenomena and in this 
sense, we have adopted the very simplified assumption that the interstellar medium is 
a constant density medium with riism = '^/cm^. Consequently, the afterglow emission 
obtained is very smooth in time. We are now interested in seeing if in this frameowrk we 
can also explain most of the time variability observed by BATSE, all of which, except 
for the P-GRB, should correspond to the beam-target phase in the above paradigm. 

We first recall the constitutive equations (Ruffini, et al. (1999,2000) [19, 20], Bianco, 
et al. (2001) [2]): 



/ AM- r2 /AM- (-2^ ^ 
Ai^int = Pb^Vi\ 1 + 271 h 



(3 a) -pBiVi 1 



71 + 



1 + 271^ 



P-Bl Vi 



(4) .; = (i+.)(,-i) = (i + .)(«- JM'''f+-^ ) 

where the quantities with the index "1" are calculated before the collision of the ABM- 
pulse with an elementary ISM shell of thickness Ar and the quantities with "2" after 
the collision. We indicate by A_Eint the increase in the proper internal energy due to 
the collision with a single shell and by ps the proper energy density of the swept up 
baryonic matter, by V the ABM pulse volume in the comoving frame, by Migm the ISM 
mass swept \vp until radius r in laboratory frame and by 7 the Lorcntz factor of the 
expanding ABM pulse, is the arrival time of the signals on the detector, counting 
from the arrival of the first photon, z is the cosmological redshift of the source and t is 
the emission time of the signal, counting from the dyadosphere formation. Details are 
given in RuHini, et al. (2001f) [24]. 

In order to proceed, we first distinguish two different regimes in the afterglow (see 
Fig. 2): in the first the intensity of the afterglow increases with time, in the second 
it decreases. The first regime goes from point 4, corresponding to the emission of the 
P-GRB (see Fig. 1), to the point P, where the peak of the radiation of the afterglow 
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is reached. During this regime, the amount of material engulfed from the interstellar 
medium is too small compared to the initial kinetic energy of the ABM pulse and the 
Lorcntz gamma factor is slowly decreasing with time, so much so that wc can assiimc 7 is 
constant in this regime. The flux emitted by the afterglow is given (RufEni, et al. (2001f) 
[24]), as a function of the laboratory time t by 



(5) F OC JO^Uismt'^, 

where 70 is the value of the Lorentz gamma factor at the moment of transparency. This 
expression can be simply expressed in terms of the arrival time ta (Ruffini, et al. (2001f) 
[24]) 



(6) F (X JO^Uism {tt 



The second regime occurs as soon as the mass-energy accreted from the interstellar 
material is no longer negligible with respect to the initial kinetic energy of the ABM 
pulse (Ruffini, ct al. (2001f) [24]). The flux emitted by the afterglow decreases now with 
the laboratory time following the law 

(7) F CX JpHp%ismt-^, 

where 7p is the value of the Lorentz gamma factor at the point P and tp is the value of 

the laboratory time when the point P is reached. There arc actually two different peaks in 
the radiation flux, if the phenomenon is the spike in the laboratory frame or in the frame 
of an asymptotic observer comoving with the detector (see RufBni, et al. (2001f) [24] for 
details). Here we consider the peak in the laboratory time. The peak of the radiation 
occurs at a value of 7 given by (Ruffini, et al. (2001f) [24]): (Mj^^/Mb) ~ 10"^ and 

(8) 7P - 0.6770, 

where Mb is the initial baryonic mass of the ABM pulse and Mj^m is the mass of the 

ISM engulfed by the ABM pulse at the time tp. Again wc can express the energy flux 
given in Eq.(7) as a function of the arrival time as (Ruffini, et al. (200 If) [24]) 

(9) F<xnism{tty^'\ 

in very good agreement with the results of the BATSE, RXTE and Chandra satellites 
(see Letter 2). The corresponding diagrams are summarized in Fig. 2. 

Once the two results presented in Fig. 1 and Fig. 2 have been understood, we can 
proceed to attack the specific problem of the time variability observed by BATSE. 

The fundamental point is that in both regimes the flux observed in the arrival time 
is proportional to the interstellar matter density: any inhomogcncity in the interstellar 
medium Auism /nism will lead correspondingly to a proportional variation in the intensity 
AI/I of the afterglow, which can indeed be erroneously interpreted as a burst originating 
in the "inner engine" 

There is a very significant signature of this kind of intensity contrasts: the AI/I is 
independent of the special moment of observation during the afterglow era and is only 
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Fig. 4. - The density contrast of the ISM cloud profile introduced in order to fit the observation 
of the burst of GRB991216. The dashed line indicates the average uniform density n = lcm~^. 



function of the density contrast. In particular, for the main burst observed by BATSE 
(see Fig. 3a) we have 

(10) (A7/T) = {Ariism/nism) ~ 5. 

There are still a variety of physical circumstances which may lead to such density inho- 
mogcncitics. 

The additional crucial parameter in understanding the physical nature of such in- 
homogeneities is the time scale of the burst observed by BATSE. Such a burst lasts 
Ata — 20,s and shows substructures on a time scale of ^ l.s (see Fig. 3a). In order to 
infer the nature of the structure emitting such a burst we must express these times scales 
in the laboratory time (see Letter 1). Since we are at the peak of the GRB we have 
7p 159 (see Eq.(8)) and Ata corresponds in the laboratory time to an interval 

(11) At = 7|, X Ata ~ 7.5 X lO^s, 

which determines the characteristic size of the inhomogeneity creating the burst AL ~ 
2.2 X lO^^cm. 

It is immediately clear from Eq.(lO) and Eq.(ll) that these are the typical dimensions 
and density contrasts corresponding to a small interstellar cloud. As an explicit example 
we have shown in Fig. 4 the density contrasts and dimensions of an interstellar cloud with 



ON THE PHYSICAL PROCESSES WHICH LIE AT THE BASES OF TIME VARIABILITY OF GRBS 9 



4.5e+045 
4e+045 

3.5e+045 
3e+045 

'eo' 

3. 2.5e+045 

.i 2e+045 

E 

_J 

1.5e+045 
1e+045 
5e+044 



a) 



100000 



Detector arrival time (tg ) (days) 

4 5 6 



With ctianging njj,^ 
=1 cm"^ 



With fixed rii. 



200000 



300000 



400000 



500000 



600000 



Detector arrival time (t, ) (s) 



700000 



2.5e+044 



2e+044 



1.5e+044 



1e+044 



5e+043 



Detector arrival time (t^ ) (days) 
25 30 35 40 




1e+006 



1.5e+006 



2e+006 



2.5e+006 



3e+006 



3.5e+006 



4e+006 



4.5e+006 



Detector arrival time (t^ ) (s) 



Fig. 5. - a) Same as Fig. 3b with the ISM cloud located at a distance of 3.17 x lO^'^ cm from the 
EMBH, the time scale of the burst now extends to ~ 1.58 x lO^s. b) Same as a) with the ISM 
cloud at a distance of 4.71 x lO^^cm from the EMBH, the time scale of the burst now extends 
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an average density < n >= l/crn?. Such a cloud is located at a distance of ~ 8.7x lO^^cm 
from the EMBH, gives rise to a signal similar to the one observed by BATSE (see Fig. 3b). 

It is now interesting to sec the burst that would be emitted by the interaction of 
the ABM pulse with the same ISM cloud if it were encountered at later times during 
the evolution of the afterglow. Fig. 5a shows the structure of the burst at a distance 
2.59 X lO^^cm, corresponding to an arrival time delay of ^ 2 days, where the Lorentz 
factor is now 7* ~ 3.60. Although the overall intensity is smaller, the intensity ratio of 
the burst relative to the average emission is consistent with Eq.(lO), but the time scales 



of the burst are longer by a factor ( ^ ) — 2 x 10^. Fig. 5b shows the corresponding 



quantities for the same ISM cloud located at a distance 3.9 x 10^ cm from the EMBH, 
corresponding to an arrival time delay of ~ 1 month, where the Lorentz gamma factor 



The approximations adopted in this paper in the solution of Eqs.(3,4) have been 
explicitly presented in all details in Ruffini, et al. (2001f) [24]. 

It is then clear that all the fundamental information on relativistic astrophysics about 
the EMBH dyadosphoro as well as its formation during the process of gravitational col- 
lapse have to be inferred from the data on the propereties of the P-GRB (Bianco, et 
al. (2001) [2], RufSni, et al. (2001g) [25]). 

The data on the E-APE appear to give mainly information on the structure of the 
ISM clouds in star-forming regions in far away galaxies. 

It is then possible to carry out, very efficiently, the sort of problematic examined, 
within our own galaxy, by the BeppoSAX satellite (see Bocchino & Bykov (2000) [4] and 
references therein). In these works the interstellar clouds have been examined using as 
"the beam" the material ejected in supernova remnants, and as "the target" a variety 
of ISM clouds in our galaxy. By properly taking into account the results summarized in 
Fig. 1 and Fig. 2 it is in principle possible, using different GRBs, to map the interstellar 
matter distribution in star- forming regions in far away galaxies at arbitrary red shift. 

This leads us into the domain of another science, of classical astronomy, into which 
the object of this work does not allow us to go today. 
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